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Abstract 

Solutions to EYM systems in 5 spacetime dimensions possessing no gravity decoupling 
limits, feature a peculiar critical behaviour which is absent in their 6, 7 and 8 dimensional 
counterparts which do possess flat space limits. This critical behaviour in 5 dimensions 
persists even when a scalar matter field is added, rendering the model nontrivial in the 
gravity decoupling limit. To this end, both regular and black hole spherically symmetric 
solutions to the higher curvature EYM-Grassmannian sigma model model in d = 5 space- 
time dimensions are constructed. A study of the solutions to the Grassmannian model in 
flat space is also carried out. 
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1 Introduction 



Gravitational theories in higher dimensions are of current interest in the contexts of (10 di- 
mensional) superstring theory and of theories with large and infinite extra dimensions with 
non-factorisable metrics. It is thus interesting to study the properties of the corresponding 
Einstein- Yang-Mills (EYM) systems in higher dimensions generalising the usual four dimen- 
sional EYM model whose regular solutions were constructed in [Jj, and black-hole solutions in 

mm 

Since the low energy effective action of string theory includes higher order terms in the 
gravitational and Yang-Mills (YM) curvatures, such systems were recently studied in spacetime 
dimensions d = 6,7, 8 j3] and d = 5 in which dimensions it is necessary to include higher 
curvature YM terms to enable the existence of particle like solutions. 

As expected, the gravitating YM models in spacetime dimensions d = 5, 6, 7, 8 support 
particle like solutions for a finite range of the gravitational parameter a^, say to a'^ax- This 
is similar to the case of gravitating monopoles [HI IZI in c? = 4. What actually happens here 
differently from the latter case is that two solutions exist for a given value of a^. But there 
is a marked difference between the ci = 6, 7, 8 cases |4j, where there exist two solutions for all 
values of a^, and the d = 5 case [S] where this is true for values of up to a critical value 
a^. In the latter case, solutions exist for values of oscillating about a^. The tracking of this 
peculiar singular behaviour in the d = 5 model is the aim of the present work. 

Now the most obvious difference between the solutions of the models |1] in = 6, 7, 8 on 
the one hand, and those of the model [5] in d = 5 on the other, is that the solutions in the 
former cases persist in the gravity decoupling limit ^, while those in the d = 5 model do not ^. 
We have therefore added a scalar matter field to the d = 5 model, which renders the gravity 
decoupling limit nontrivial, to investigate the nature of the critical behaviour discovered in 
[HI. We thus answer the question: Is this critical behaviour a consequence of the absence of a 
gravity decoupling limit in c? = 5, or is it a peculiarity of the dimensionality of the spacetime 
itself? Our answer is, that this property pertains to the dimensionality of the spacetime, since 
we find that it persists also in the new model which supports a gravity decoupling solution. 

For this purpose we introduce an SU (2) gauged 4x2 Grassmannian field ^Uj describing a 
sigma model, to the d = 5 EYM model. The introduction of this scalar field is analogous to 
the inclusion of a Higgs^ field in d = 4 as in Refs. E] . 

"'^The solutions in the flat space limit of the d = 6,7,8 models 0] are exemplified by the instanton-like 
solution studied in 8,. 

^Note that due the Derrick scaling requirement the YM field in the static 4 Euclidean dimensions supports 
a 'soliton' only if the YM system consists of the p = 1 term, exclusively. On the other hand when gravity is 
switched on, the absence of the p = 2 YM term (see Ref. [0] for the YM hierarchy) prevents the existence of a 
soliton, due to the same scaling requirement. 

•^The choice of a Grassmannian field rather than a Higgs field is because the gauge connection in a d — 1 
dimensional Higgs model ^21 supporting a 'soliton' behaves as one half pure gauge asymptotically, resulting in 
decay. As a result the integral of the p = 1 YM term is convergent only in (d — 1) < 3 so that only higher p 
YM terms are admissible for (d—l) > 4. So if we insist in keeping the usual p — 1 YM term, then we must avoid 
using of a (generalised) Higgs model ^21 ■ This contrasts with the faster decay of an 'instanton' in a pure gauge 
theory where the gauge connection is pure gauge and hence the integral of the p = 1 YM term is convergent 
in c? — 1 = 4. It turns out that the connection of the gauged Grassmannian model LlUj is asymptotically pure 
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We find that the singular behaviour in question, which is pecuhar to the solutions in d = 5 
spacetime only, presists whether or not the model supports a regular solution in the flat space 
limit. 

2 The model and the equations 

In the first subsection we give the Lagrangian of our 4+1 dimensional model in Minkowski space, 
and in the second one we impose static spherically symmetry and write down the resulting one 
dimensional ordinary differential equations. 

2.1 The models 

We take the gravitational and YM sectors of the model to be precisely the one analysed in [S] , 
augmented by the Grassmannian sigma model term 

^ ^grav ~\~ ^YM ~\~ ^ grass ; (1) 

in 5 spacetime dimensions. The first two terms of (0) describe the EYM sector defined as 

Cgrav = e ^ , CYM = e Tr F{2f + ^ Tr ^(4)^) , (2) 

where R{i) describes the usual Einstein gravity, and e = det e° = ^J—det g^^, ej^ are the 5-beins. 
F{2p), for p = 1,2, is the 2p-form YM curvature (see Refs. jl| and which for p = 1 is the 
usual 2-form YM curvature taking values in the antihermtian representation of the algebra of 
SU{2) here. 

The Grassmannian sigma model part of (Q) is 

Cgrass = e Tr (^tD^.z^D'^z + ^ V(I + z^Tr^z) + A(I - z^ z)^ (3) 

described by the 4x2 Grassmannian field subject to the 2x2 condition 

z^z = 1 

and whose SU (2) covariant derivative is defined by 

D^z = d^z-zA^ . (4) 

In the constant /i has the dimensions of mass and leads to the exponential localisation of 
the ensuing topologically stable lump, and the Lagrange multiplier A is a 2 x 2 array. The 
gravity decoupled version of this model is the gauged Grassmannian sigma model for which the 
'instanton' solutions were constructed numerically in ^|, which will be studied in more detail 
here. 

gauge, and hence our choice. 
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2.2 The classical equations 

In d dimensional spacetime, we restrict to static fields that are spherically symmetric in the 
d — 1 spacelike dimensions with the metric Ansatz 



ds^ = -a{rfN{r)dt^ + NirY^dr"^ + r'^d^l^ 



(5) 



where r is the spacelike radial coordinate and dVLfi-2 is the d — 2 dimensional angular volume 
element. 

We take the static spherically symmetric SU (2) YM field in 5 spacetime (i.e. 4 Euclidean) 
dimensions, in one or other chiral representation of S'0±(4), to be 



An = , Ai 



1 — w 



1 /1±75 



(6) 



where 



1, 



4 [* Jl 



where in a more familiar notation S-^'* = Ui = I) and S^-"-' = 5"j = (— ia, I), with i = 
1, 2, 3, 4 and in terms of the three Pauli spin matrices a. 

The spherically symmetric Ansatz for the Grassmannian field z, whose consistency has been 
checked, is 



sin I I 

/ 

2 



(7) 



Subjecting Q and Q to spherical symmetry by empoloying the Ansatze (0), (jHI) and (j?!), and 
subjecting the resulting one dimensional Lagrange density to the variational principle, we find 
the following equations for the functions /(r), w{r), N{r) and cr{r), 

{r^aNfy + ra {3w - fih-^)smf = 0, (8) 



Ti {{raNw')' -2r-^(r{w^ - l)w) + 3r2(w;^ - 1) {r'^aN{w^ - 1)10')' = 2r r + cos /), (9) 



m 



I ^1 



Nw" + 



w 



Nw 



+ — r- 
12 



A^/'2 + — + 2m; cos / + l) + 2y?{\ - cos/) 



(10) 



The ADM mass M = lim^^oo ''^('")5 with m(r) defined as 

m(r) = n^^Kir^{l - N) . (12) 

For N{r) = a{r) = 1 everywhere, namely the gravity decouphng hmit, (jH]) and satisfy the 
flat space solution. 

In the numerical work below, we set the self interaction potential of the Grassmannian 
field, = without changing the qualitative nature of the solutions. This is because in this 
model, as in Higgs models, the finite energy condition leads to a unique asymptotic value for 
the matter field, as can be seen easily by inspection of the term multiplying r in (fTn|) . This 
situation contrasts with that in certain other gauged nonlinear sigma models, where a unique 
asymptotic value for the matter field can be ensured only by the inclusion of a self interaction 
potential, e.g. the pion mass potential of the Skyrme model, leading to bifurcations 



2.3 Boundary values and asymptotic behaviour 

In the next section, we will solve the above equations with the appropriate boundary conditions 
for the radial functions m(r), a{r) w{r) and /(r) which guarantee the solution to be regular 
at the origin and to have finite energy. For regular solutions, the boundary conditions at the 
origin are 

m(0) = , wiO) = l, /(0) = 7r, (13) 
while the conditions satisfied on the event horizon r = rh are 

N{rh) = 0, cr{rh) = an, w{rh) = Wh, f{rh) = fh, (14) 
with ahyWh, fh real constants. The asymptotic form of the solution is 

lim cr(r) = 1 , lim w{r) = —1 , lim /(r) = 0. (15) 

The condition on a{r) results in the metric being asymptotically Minkowskian. 

The asymptotic solutions to these functions can be systematically constructed in both re- 
gions, near the origin (or event horizon) and for r ^ 1. Defining = ^ we find for r ^ 1 

/(r) = TT - Car + o(r^) , 

w{r) = 1 + Cir^ + o(r^) , 

a{r) =ao[l + 2a^cy{n + 12T2cl]+o{r^)) , 

m(r) = V[c?(ri + 6c2r2) + ^rc^] +o(r6). (16) 

For black hole configurations, the expression of the solutions near the event horizon is 
fir) = /, + f(r,)(r-r,) + 0((r-r,)2), 
w{r) = Wh + w'{rh){r - rh) + 0{{r - r-hf), 

<y{r) = a^a' [l + ((n + ^^' )w'\rh) + IrrJ'^irh)) (r - r^)) + 0((r - r,)^). 



r? 



^^^) = + ^'irh){r -Th) + 0{{r -rt)^), (17) 
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where 



rrt 



m'irh) = ^-^ + ( -{wl + 2wh cos fh + I) + 2fi\l- cos fh 

12 \rl 



w'{rh) 



2rrh{wh + cos/fe) + 2TiWh{wh -l)/rh 
{riTH + St2{wI - l)yrl){-8a^m'{rj,)/rl + 2^) 



rli-8a^m'ir,)/rl + 2/r,)- ^ 



For r ^ 1 we find 



ria2c|(c?-3)^ 



a(r) 



{2d - 4)r 



2d-4 



= "^°o - g^,-l (19) 

In the first member of ()19|). K2{fir) is a Bessel function leading to exponential decay by virtue 
of the mass term in the potential. The constants ci,C2,C3,ao, fh,Wh,(Jh and moo have to 
be determined numerically; moo is nothing else but the ADM mass of the solution as noted 
previously. They depend generically on the coupling constants of the theory. 



3 Numerical results 

In this section we will present three sets of results. These are the non gravitating SU{2) gauged 
Grassmannian solitons, the gravitating regular solutions of this system, and, the corresponding 
black hole solutions. Each will be reported in a separate subsection. 

The numerical integrations were carried out using both a shooting method as well as ap- 
plying the numerical program COLSYS ^1], with complete agreement to very high accuracy. 

3.1 SU{2) gauged Grassmannian solitons 

Since an important reason for considering the model in this work is that it supports a non 
gravitating fiat space limit, it is worth studying the properties of the latter for its own sake. 
This is especially so since we find that this (non gravitating) model supports solutions which 
exhibit nodes in the profile of the Grassmann function /(r) defined in (|7j). 
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Figure 1. The functions w{r^ and /(r) and the energy density e(r) are plotted as functions of 
radius for typical flat space solutions with the coupling constants t\ = 3t2 = 1,t = 0.1. The node 
number k of the Grassmanian function /(r) is also indicated. 

The profiles of the YM functions wir) do not change appreciably for the solutions with differ- 
ent number of nodes of /(r), neither qualitatively nor quantitatively, exhibiting only one node. 
Also, somewhat unexpected, the total mass/energy of these solutions stays almost constant, 
when increasing the node number of /(r) (with differences less than one percent). We shall see 
below that all the multinode solutions result in qualitatively very similar gravitating solutions, 
which we will exploit to simplify the numerical work section ()3.3|) . 

It turns out that there are solutions with infinitely many nodes in the function /(r), such 
that the profiles contract towards the origin as the number of nodes n increases. The profiles 
of w[r) and the energy density remain insensitive to the increase in the number of nodes n of 
/(r), as seen in Figure 1. In the limit n ^ oo the profile of /(r) foo{f) shrinks to the origin 
such that 

foo{r) = for all r . 

We refer to this configuration as that of frozen /, namely / = everywhere. Also, for all 
considered solutions, as well as the gravitating counterparts, we find that the gauge function 
w{r) monotonically decreases towards its asymptotic value without presenting local extrema. 

3.2 Regular gravitating solutions 

We numerically integrate the Eqs. (jHj)- (fTT|) with the boundary conditions (fTS|) . (fT3j) for ti = 
1,T2 = 1/3 and several values of r, finding the following picture. First, for being small 
enough, a branch of solutions smoothly emerges from the fiat space configurations. When 
increases, the mass parameter M decreases, as well as the value cr(0) and the minimum Nm 
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Figure 2. The value Nm of the minimum of the metric function N , the value of the metric function 
a at the origin ct(0), as well as the mass M are shown as a function of = nd/{8Ki) for regular 
solutions with ri = 3t2 = 1 and two different values of r. The mass of the flat space solution is 
M = 1.00236 for r = 1 and M = 0.49257 (r = 0.1). 

of the function N{r) decrease, as indicated in Fig. 2. These solutions exist up to a max- 
imal value amax of the parameter a, which is smaller than the corresponding value in the 
pure EYM theory 0, and depends on the value of the coupling parameter r. For example, 
we find numerically af^ax ~ 0.2573 fror r = 0.1 while the corresponding value for r = 1 is 
'^max ~ 0.06855. (Without a Grassmanian field, this branch extends up to afnax ~ 0.5648.) 

Similar to the EYM case we found another branch of solutions on the interval G 
['-'^cr(i)' '^max] with depending again on the value of r (e.g. a^r{i) ~ 0.06302 for r = 1). 

On this second branch of solutions, both cr(0) and Nm continue to decrease but stay finite. 
However, a third branch of solutions exists for G [a^^^^), ^^(2)] ' which the two quantities 
decrease further. A fourth branch of solutions has also been found, with a corresponding a^r(3) 
close to C('^r(2)- Further branches of solutions, exhibiting more oscillations very likely exist but 
their study is a difficult numerical problem. Progressing on this succession of branches, the 
main observation is that the value a{0) decreases much faster than that of Nm as illustrated 
in Fig. 2. The pattern strongly suggests that after a finite (or more likely infinite) number of 
oscillaions of cr(0), the solution terminates into a singular solution with cr(0) =0 and a finite 
value of N{0). As seen in Figure 2, the mass parameters do not increase significantly along 
these secondary branches. 

This is the behaviour observed in [5| for the pure EYM theory. The inclusion of a Grass- 
manian extra field does not seem to qualitatively change the properties of the system. 

In Fig. 3, we present the profiles of the metric functions and cr for the same value of 
on the first and second branch. 
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Figure 3. The metric functions N(r),a{r) are plotted as functions of radius for typical regular 
gravitating solutions. 

The discussion in this section is restricted to the Grassmanian field presenting no nodes. We 
do not expect the consideration of solutions with nodes to change this picture. 



3.3 Black hole solutions 

According to the standard arguments, one can expect black hole generalisations of the regular 
configurations to exist at least for small values of the horizon radius Vh- This is confirmed by 
the numerical analysis for solutions with no nodes in /(r) as well as solutions with a frozen 
Grassmanian field /(r) = foe = 0. 

Again, the properties of the solutions we find are rather similar to the five dimensional black 
hole solutions without a Grassmanian field discussed in [S]. Firstly, black hole solutions seem 
to exist for all values of a for which regular solutions were constructed. Also, for a given set of 
couplings (ri,r2,r), the solutions exist only for a limited region of the {rh,a) space. 

The typical behaviour of solutions as function of Vh is presented in Figure 4, for a small 
value of a as compared to to maximal value amax of the regular solutions. Starting from 
a regular solution and increasing the event horizon radius, we find a first branch of solutions 
which extends to a maximal value r^rnax)- The variation of mass and cr(rft) is relatively small on 
this branch. Extending backwards in rh, we find a second branch of solutions for < rh{rnax)- 
This second branch stops at some critical value r^cr), where the numerical iteration fails to 
converge. The value of a{0) on this branch decreases drastically, as shown in Fig. 4. Also, the 
surface gravity k of the solutions, given by 
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Figure 4. The value Wh of the gauge function at the horizon, the value ah of the metric function a 
at the horizon, the mass M, as well as the surface gravity k divided by 12 are shown as functions of 
the event horizon radius for black hole solutions with = 0.1, ri = 3t2 = 1, r = 0.1. These results 
are obtained for a "frozen" Grassmanian field / = 0. 

strongly decreases on this branch, approaching a very small value. However, the increase 
of the total mass is still very small. Similar to the EYM case jS], higher branches of solutions 
on which the value o"(0) continues to decrease further to zero are likely to exist. However, the 
extension of these branches in will be very small, which makes their study difficult. 

Although the results in Figure 4 correspond to a frozen Grassmanian field, we do not expect 
a different result for solutions with a nontrivial f{r). Similar to the regular case, we find that 
the node number of the Grassmanian field does not significantly affect the propertiess of the 
gravitating solutions. 

However, the global picture we find (and the corresponding EYM results) may change by 
considering values of a near amax- We hope to come back on this point in a future work. 

4 Summary and discussion 

The aim of this work is to find out whether a particular singular behaviour of solutions of a 
gravitating YM model^ in = 5 spacetime is a persistent feature of the dimensionality of the 
spacetime? The reason for asking this question is that the EYM model in question here does 
not feature this critical behaviour in spacetimes ci = 6, 7, 8. 

More simply stated, the model(s) in c/ = 6, 7, 8 exhibit only a maximum value of the grav- 

^The definition of the model in question is d (spacetime) dependent. While in all d these are formally the 
same, nevertheless the gauge groups are different, depending on d 
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itational coupling = a^^^, while that in c? = 5 has, in addition to an = af^^^, also an 
a^j. > 0. This is seen from the oscillatory behaviour of converging to a^^ in Figure 2. Because 
d = 5 spacetime is particularly relevant in the context of the AdS/CFT correspondence, and 
because the solutions in c? = 5 spacetime differ markedly in this respect from those in c? = 6, 7, 8 
spacetimes, it is important to study this peculiar feature in d = 5 further. 

The original model(s) introduced in j3] involve a dimensionful constant in addition to the 
gravitational constant, analogously with the d = 4 gravitating YMH model [HI HI- Like the 
latter, regular solutions to the higher dimensional EYM models exist only for values of the 
gravitational coupling up to a maximum a^^^. But this analogy with the EYMH case goes 
further only in the d = 5 EYM case, where in addition to a"^^^, there occurs also a critical 
a'^j. > 0. In the d = 6, 7, 8 EYM cases it appears that = 0. This is surprising since the 
d = 6,7,8 EYM models support regular solutions in the flat space limit just like the d = 4 
EYMH model, while the d = 5 EYM model does not support a fiat space solution. 

It is to throw some light on this question that we modified the d = 5 EYM model by 
introducing a (Grassmannian) scalar matter field, which results in the new model supporting 
a fiat space solution. The result is that the qualitative features of the solutions of the original 
d = 5 EYM model are preserved. 

In passing, we studied the (static) solitons of the d = 5 fiat space SU (2) gauged Grassman- 
nian model, and found that these form an infinite sequence of solutions exhibiting multi-nodes 
in the profile of the Grassmannian profile function. 
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